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INTRODUCTION 


Nonlinearities now appear regularly in all phases of aerodynamic 
research. Experimenters in dynamic stability have particular reason 
to be concerned about them, since the methods available for extract- 
ing stability derivatives from dynamical data for the most part do 
not acknowledge the presence of nonlinearities. When the aerodynamic 
characteristics are nonlinear, these methods continue to give 
results, but it becomes very difficult to say precisely what the 
results mean. Consider, for example, the two widely used experi- 
mental methods sketched below. 




In the first method, the sting is brought to a fixed mean angle of 
attack ccjn and the model is forced to oscillate harmonically about 
with a small amplitude a 0 . The work done per cycle is meas- 
ured and identified with the aerodynamic damping derivative 
Cmq + In the second method, the model is displaced from zero 

angle of attack against the resistance of a spring and then released. 
The decay of the free oscillation about zero is measured and iden- 
tified with the damping derivative Cm^ + Cm^. When the aerodynamic 
moment is a linear function of its variables, the same value for the 
damping derivative is obtained from both experiments. When the 
moment is nonlinear, two entirely different results can be obtained. 
Neither result is necessarily wrong; each merely reflects the par- 
ticular method in use, which in turn reflects a particular facet of 
the underlying nonlinear phenomenon. 

Consider the analyst who must use these results to analyze 
the motion of a vehicle in flight. If the motion he anticipates 
consists of small deviations from a trim condition, the first result 
may be applicable. If the motion he anticipates consists of large- 
amplitude oscillations around zero, the second result may be 
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applicable. But, suppose he anticipates the second and has results 
only for the first. Or, suppose he can anticipate neither motion. 

He has no assurance at all that the results at hand have any mean- 
ing in his problem. Obviously, this situation imposes a severe 
limitation on the usefulness of wind-tunnel measurements of stability 
coefficients . 

It is clear that if this limitation is to be overcome, results 
from different facilities and different experimental techniques must 
be made to correlate , so that they may be used interchangeably. In 
other words, the experimental results must be freed from dependence 
on the particular method used to determine them. A recent report 
(ref. 1) addressed to this problem is briefly summarized in the 
remainder of this note. 


GENERAL FORM FOR THE AERODYNAMIC PITCHING MOMENT 


Before the results of an experiment can be freed from the 
method used to determine them, the true form of the equation govern- 
ing the motion must be established. This form will, of course, 
depend on the form of the aerodynamic forces and moments. The prob- 
lem is, therefore, to write the aerodynamic forces and moments in a 
form which is sufficiently general to apply to any of a wide variety 
of possible motions. In the linear case, one achieves this by writ- 
ing the forces and moments as a sum of terms involving stability 
derivatives. The question then arises, do stability derivatives 
even exist in the nonlinear case? Based upon the analysis of ref- 
erence 1 the answer is yes , with certain strict reservations . For 
the longitudinal case, involving arbitrary variations of angle of 
attack a and angle of pitch 9 , the aerodynamic pitching moment 
has the following form 

C m (t) = C m (co;a(t) ) + £(t) —■ Cm^ocqaCt)) 

+ a(t)“ C ffl(i (a(t)) (1) 

The terms have the following meaning: The first term C m («>;a(t)) 

is the familiar steady -state pitching -moment coefficient due to 
angle of attack, evaluated at each instant as though the instan- 
taneous value of a, a(t) , were fixed for an infinite time at that 
value. The infinity symbol as used in this notation is merely a 
reminder that the flow is fixed for all time, that is, is steady 
at the particular value of a under consideration. The second 
term Cm^(oo;a(t) ) is the rate of change with 9 of the pitching- 

moment coefficient that would be measured in a steady flow , 
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evaluated at 0=0 with a fixed for an infinite time (thus the 
infinity sign again) at the instantaneous value a(t) . The third 
term Cm& is defined by means of an integral: It is the area 

bounded by the indie ial pitching -moment curve and the final value 
of the curve, the indicial response being evaluated with a fixed 
at the instantaneous value a(t) and with § fixed at zero. 

It must be emphasized that the definitions just given for the 
terms in equation (l) are precise and unique; they arise from a 
rigorous analysis correct to the first order in frequency. Within 
this order of approximation, alternative forms or definitions would 
not be rigorously justifiable. Thus, for example, note that each 
of the terms in equation (l) shows a dependence only on the instan- 
taneous angle of attack a(t); to show a dependence on a(t) , 5(t), 
or higher derivatives of a and 0 would not be justifiable. 
Equation (l) is the desired general form which should underlie all 
special cases and apply to all particular motions. 


FORM OF THE AERODYNAMIC PITCHING MOMENT IN A SPECIAL CASE 


Further specifications of the form of C m (t) are possible in 
particular cases. Let us consider one of these for a wind-tunnel 
experiment. In most wind-tunnel experiments, the model is pinned 
to a fixed point at its axis of rotation. Then 5 = &. For the 
wind-tunnel experiment, therefore, Cmq and Cm^ simply add, and 
Cm(t) takes the form 

C m (t) = c m( 00 J <x ( t )) + ^(t) [CmqO^cc-Ct)) + C^aCt) ) ] (2) 

Now suppose that the model under study is known from static tests 
to have a static pitching -moment curve which is representable as 
an odd cubic in a over a substantial range of a. Then C m (oo;a(t)) 
has the form of equation (3) below, and the rate of change of C m 
with a, Cnu(°°;a(t) ) , has the form of an even quadratic in a 
(eq. (10). 

Cm(co;a(t)) = a(a + ba 2 ) (3) 

C^CoojaCt)) = a + 3b<x 2 (4) 

Since the boundary conditions which yield C^ are very similar to 
those which yield Cr^, it is very reasonable to assume that C^ 
likewise will show no more than an even quadratic dependence on a. 
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under similar conditions, it can be argued that 0^ also will be 
an even quadratic in a. Then, so is the sum + Cm^- Thus, 

Cm^C 00 ; 0 ^)) = c + da 2 (5a) 

Cny^aCt)) = e + fa 2 (5b) 

and 

C ®q + Cm a = h o + h 2 a2 (5c) 

Alternatively, can be written as a linear function of the 

steady -state parameter C n i a (cc;a(t) ) as shown below. 

C'mq "b = C + DCm <x (o o ;a(t) ) (6) 

This particuleirly simple form for the damping coefficient has certain 
consequences which could be checked experimentally. 

EXPERIMENTAL CORRELATIONS 


Now let us consider the main objective, correlation of 
experimental results. We wish to show that, at least in certain 
cases, knowledge of the form of the pitching -moment coefficient 
enables one to free the results for damping coefficient obtained by 
a particular experimental method from dependence on that method. 

If the results from two very different experimental methods were 
both freed from dependence on their respective methods, then those 
results should correlate. Results which correlate can of course 
be used interchangeably. If results from two very different motions 
can be used interchangeably, then we are assured that either result 
is a general one; that is, one which is applicable to any motion 
whatever. 

Let us assume that the model under consideration fulfills the 
requirement discussed in the last section, namely, that its static 
pitching -moment curve be a cubic in a. Then we anticipate that 
the damping coefficient should be a quadratic in a or, alterna- 
tively, a linear function of the static pitching -moment curve 
slope Cn^oo^t) ) • Suppose this is true. Then this is the form 
we expect whenever a particular experimental result for damping 
coefficient is freed from the method of obtaining it. We next 
consider how this might be achieved in specific experiments. 
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EXPERIMENIAL TECHNIQUES 


Consider the two very different experimental methods illustrated 
in the previous sketch. The first, forced oscillations about a fixed 
mean angle of attack <% , is shown again below. 


a a 0 




The work required to drive the model is measured and equated to the 
work done by the aerodynamic damping moment. In the linear case, 
the work measured can be identified with a single, unique value of 
the aerodynamic damping coefficient Cm^ + C^. In the nonlinear 

case, this is no longer true. Here, when the work measured is 
identified with the work done by aerodynamic damping, the result 
obtained depends on both the mean angle and the oscillation 

amplitude a 0 ; that is, the result depends on the method used to 
obtain it. The objective then is to free the result from this 
dependence. 

When the aerodynamic damping coefficient depends linearly on 
the instantaneous value of the static pitching -moment curve slope 
(i.e., eq. (6)), the work done by damping will have the form shown 
in equation (7) where C and D are the same C and D in equation 
(6) . 


where 


Work/Cycle = C + DC mae 



( 7 ) 


(°°;<Xm) 


( 8 ) 
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That is, the work done by damping becomes a linear function of an 
effective value of 0^, where Cn^ can be evaluated from 

static data and equation (8) . Hence, plotting the work done by 
damping against gives a straight line. 



Note that the same straight line will be obtained for any value 
of oscillation amplitude a Q . The slope and zero intercept of the 
line give the values of C and D. When these are substituted in 
equation (6), together with the true static curve Cm a (oo;a(t) ) , the 
form obtained for the damping coefficient C m ^ + should then be 

general - applicable to any motion whatever. This expectation can 
be checked by an independent experiment . 

The second method, to be used as a check on the generality of 
the preceding result, is the well-known free-oscillation technique, 
shown again below. 
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The equation of motion for this system is 


la 

<lo Si 


= C m (t) + spring terms 


(9) 


With the same model as the one used in the first experiment , the 
static pitching -moment curve is a cubic in a and the damping 
coefficient is a quadratic in a, so that the pitching moment 
(eq. (2)) takes the form 

C m (t) = a(a + ba 2 ) + — (h Q + hsa 2 ) (10) 

v o 

Combining equations ( 9 ) and (10) yields the following for the equa- 
tion of motion 



+ ba 2 ) 


+ 


yp- (h Q + h;ga 2 ) + spring terms 
v o 


( 11 ) 


In this form, the equation of motion is a combination of the Duffing 
equation and the Van der Pol equation. An approximate solution to 
the equation is obtainable from the Kryloff -Bogoliuboff method. Fit- 
ting the solution to the measured motion history enables the evalua- 
tion of the four constants a, b, h Q , h 2 . These have a known 
correspondence with the constants C and D used in the first 
experiment ; that is, C and D can be expressed as 



( 12 ) 


If the values of C and D computed from equations (12) agree with 
those obtained in the first experiment, then we have shown that 
two very different experimental methods can be made to yield results 
which correlate. This in turn assures that the form which corre- 
lates is a general one, applicable to any motion whatever. 


SUGGESTED EXPERIMENT 


The above discussion suggests an obvious experiment: Select a 

model whose static pitching moment C m (oo;a) can be represented as a 
cubic over a substantial range of a. Carry out the two experiments 
in the manner just described. If the results for Cn^ + Cm^ can be 


expressed as a quadratic in a, and if the values of C and D agree 
between the two experiments, this will count both as a confirmation 
of the theory and a useful means of widening v the applicability of 
the experimental results. Finally, we point out that a linear rela- 
tionship between Cm^ + and can be expected to obtain 

under a number of circumstances other than those discussed here 
(i.e., Cnu a quadratic function of a). For example, a linear 
relationship between Cm^ + and can be expected to obtain 

whatever the variation of C^, provided that Cm^ is either small 

or essentially independent of a. The sketch below shows an experi- 
mental result which confirms this statement for a model whose value 
of Cniq is known to be small. Note that is a complicated 

function of a; nevertheless, the relationship between work done and 
Cma e is quite obviously linear. 
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